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ABSTRACT 

We describe an overall picture of galactic-scale star formation. Recent high-resolution magneto-hydrodynamical simulations of two- 
fluid dynamics with cooling/heating and thermal conduction have shown that the formation of molecular clouds requires multiple 
episodes of supersonic compression. This finding enables us to create a scenario in which molecular clouds form in interacting shells 
or bubbles on a galactic scale. Eirst we estimate the ensemble-averaged growth rate of molecular clouds over a timescale larger than a 
million years. Next we perform radiation hydrodynamics simulations to evaluate the destruction rate of magnetized molecular clouds 
by the stellar EUV radiation. We also investigate the resultant star formation efficiency within a cloud which amounts to a low value 
(a few percent) if we adopt the power-law exponent ~ -2.5 for the mass distribution of stars in the cloud. We finally describe the 
time evolution of the mass function of molecular clouds over a long timescale (>lMyr) and discuss the steady state exponent of the 
power-law slope in various environments. 

Key words, galaxies: ISM - ISM: clouds - ISM: magnetic fields - ISM: molecules - stars: formation - ISM: kinematics and 
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1. Introduction 


Recent observational studies of nearby star-forming regions 
with the Herschel Space Observatory have convincingly 
shown that stars are born in self-gravitatin g filaments (e.g., 
Andre et aP 120101 : lArzoumanian et akl 1201 lb . In addition, the 
resultant mass function of star-forming dense cores are now 
expla i ned by the mass distribution along filaments (llnutsukal 
20011: lAndr6 et’aP l2014l) . This simplifies the question of the 
initial conditions of star formation, but poses the question 
of how such filamentary molecular clouds are created in the 
interstellar medium (ISM) prior to the star formation process. 
Recent high-resolution magneto-hydrodynamical simulations 
of two-flui d dynamics with cooling/heati ng and thermal con¬ 
duction by llnoue & Inutsukal (l2008lT2()0^ have shown that the 
formation of molecular clouds re quires multiple episodes of 
super s onic compression (see a lso iHeitsch, Stone. & HartmannI 
|2009|) . llnoue & Inutsukal j2012h further investigated the forma¬ 
tion of molecular clouds in the magnetized ISM and revealed 
the formation of a magnetized molecular cloud by the accre¬ 
tion of HI clouds created through thermal instability. Since 
the mean density of the initial multi-phase HI medium is 
an order of magnitude larger than the typical warm neutral 
medium (WNM) density, this formation timescale is shorter 
than that of molecular cloud formation solely through the 


accumulation of diffuse WNM (s ee, e.g., Kovama & Inutsuk^ 


2002 : Hennebelle et al. 20081: Hdtsch & HartmannI 20081 : 


Baneriee et alJ 120091 : IVazauez-Semadeni et alJ 1201IL for the 


cases of WNM flows). The resulting timescale of molecular 
cloud formation of >10 Myrs is consistent with the evolutionary 
times cale of molecular clouds in the EMC dKawamura et al.l 
120091) . 


We have done numerical simulations of additional com¬ 
pression of already-formed but low-mass molecular clouds, and 
found interesting features associated with realistic evolution. 
Figure 1 shows a snapshot of the face-on view of the layer 
created by compressing a non-uniform molecular cloud with a 
shock wave propagating at 10 km/s. The direction of the shock 
compression is perpendicular to the layer. The magnetic held 
lines are mainly in the dense layer of compressed gas. The 
strength of the initial magnetic held prior to the shock com¬ 
pression is 20//Gauss and that of the dense region created after 
compression is about 200juGauss on average. Many dense fila¬ 
ments are created with axes perpendicular to the mean magnetic 
held lines. We can also see many faint filamentary structures that 
mimic “striations” observed in the Taurus Dark Cloud and are al¬ 
most parallel to the mean magnetic field lines dGoldsmith et al.l 
l2008l) . In our simulations, these faint filaments appear to be feed¬ 
ing gas onto dense filaments (simi lar to w hat is observed f or lo¬ 
cal clouds by ISugitani et al.l (e.g., l201 ll) : IPalmeirim et al.l (e.g.. 
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Fig. 1. Face-on column density view of a shock-compressed dense layer 
of molecular clouds. We set up low-mass molecular clouds by the com¬ 
pression of two-phase HI clouds. This snapshot shows the result of an 
additional compression of low-mass molecular clouds by a shock wave 
propagating at 10 km/s. The magnetic field lines are mainly in a dense 
sheet of a compressed gas. The color scale for column density (in cm“^) 
is shown on top. The mean magnetic field is in the plane of the layer 
and its direction is shown by white bars. Note the formation of dense 
magnetized filaments whose axes are almost perpendicular to the mean 
magnetic field. Fainter “striation”-like filaments can also be seen, that 
are almost perpendicular to the dense filaments. 


l2013h : iKirk et al] (e.g., l2013h T Once the line-mass of a dense 
filament exceeds the critical value (2C?/G) star formation is 
expec ted to start (llnutsuka & Mivam^l 19921 1997t lAndre et~aD 


I 2 OIOI) . This threshold of line-mass for star formation is equiv¬ 
alent to the threshold of the column density of molecular gas 
116Mopc“^ (iLada et al.l 201 Oh. if the widths of filaments a re all 
close to 0.Ipc (lArzoumanian et al.ll201 ItlAndre et alJl2014l) . 

Although further analysis is required for quantitative com¬ 
parison between the results of simulation and observed struc¬ 
tures, Figure 1 clearly shows that the structures created by mul¬ 
tiple shock wave passages do match the characteristic structures 
observed in filamentary molecular clouds. This motivates us to 
describe a basic scenario of molecular cloud formation. The 
present paper is focused on the implications of this identifica¬ 
tion of the mechanism of molecular cloud formation. 


2. A Scenario of Cloud Formation Driven by 
Expanding Bubbles 

HI observations of our Galaxy re veal many shell-like struc¬ 
tures near the ga lactic plane (e.g., Hartmann & Burt^ 119971 : 
ITavlor et al.ll2()03i) . We identify repeated interactions of expand¬ 
ing shock waves as a basic mechanism of molecular cloud for¬ 
mation, and we depict the overall scenario of cloud formation 
in our Galaxy as a schematic picture in Figure |2] In this pic¬ 
ture, red circles correspond to the remnants of shock waves due 
to old and slow supernova remnants or expanding HII regions. 
Cold HI clouds embedded in WNM are almost ubiquitously 



Fig. 2. A schematic picture of sequential formation of molecular 
clouds by multiple compressions by overlapping dense shells driven 
by expanding bubbles. The thick red circles correspond to magnetized 
dense multi-phase ISM where cold turbulent HI clouds are embedded 
in WNM. Molecular clouds can be formed only in the limited regions 
where the compressional direction is almost parallel to the local mean 
magnetic field lines, or in regions experiencing an excessive number of 
compressions. An additional compression of a molecular cloud tends 
to create multiple filamentary molecular clouds. Once the line-mass of 
a filament exceeds the critical value, even in a less massive molecular 
cloud, star formation starts. In general, star formation in a cloud acceler¬ 
ates with the growth in total mass of the cloud. Giant molecular clouds 
collide with one another at limited frequency. This produces very unsta¬ 
ble molecular gas and may trigger very active star formation. 


found in the shells of these remnants (e.g., H artmann & BurtonI 
1 19971: iTavlor et ^ 120031: iHosokawa & Inutsuka 2007 ). Molec¬ 
ular clouds are expected to be formed in limited regions where 
the mean magnetic field is parallel to the direction of shock wave 
propagation, or in regions where an excessive number of shock 
wave sweepings are experienced. Therefore, molecular clouds 
can be found only in limited regions in shells. Note that the typ¬ 
ical timescale of each shock wave is of order IMyr, but the for¬ 
mation of molecular clouds requires many Myrs. Some bubbles 
become invisible as a supernova remnant or an HII region many 
million years after their birth. Therefore, this schematic picture 
corresponds to a “very long exposure snapshot” of the real struc¬ 
ture of the ISM. Each molecular cloud may have random ve¬ 
locity depending on the location in the most recent bubble that 
interacts with the cloud. Interestingly, this multi-generation pic¬ 
ture of the evolution of molecular clouds seems to agree with the 
observational findings of iDawson et al.l(l201 lalbLl2015l) . who in¬ 
vestigated the transition of atomic gas to molec ular gas in the 
wall o f Galactic supershells. In the case of EMC. iDawson et al.l 
(1201 3h concluded that only 12 ~ 25% of the molecular mass can 
be apparently attributed to the formation due to presently visi¬ 
ble shell activity. This may not be inconsistent with our scenario 
since Dawson et al (2013) only considered HI supergiant shells, 
whereas molecular clouds in our model can form at the interface 
of much smaller bubbles and shells (which observationally are 
more difficult to identify and characterize in the HI data) and the 
timescale for cloud forming shells to become invisible is much 
shorter than the growth timescale of molecular mass. 

A typical velocity of the shock wave due to an expanding 
ionization-dissociation front is lOkm/s, as shown by Hosokawa 
& Inutsuka (2006a), since it is essentially determined by the 
sound speed of ionized gas (~ lO'^K). Iwasaki et al. (2011b) 
have shown that if a molecular cloud is swept-up by shock wave 
of lOkm/s, it moves with a velocity slightly less than the shock 
speed. Thus, the mean velocity of each molecular cloud should 
be somewhat smaller than that of the most recent shock wave. 
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When the shock velocity of a supernova remnant is much higher 
than lOkm/s, the resultant interaction would result in the destruc¬ 
tion of molecular clouds. Therefore, the cloud-to-cloud velocity 
dispersion of molecular clouds should be similar to lOkm/s. Ac¬ 
cording to this acquisition mechanism of random velocity, the 
velocity of a cloud is not expected to depend strongly on its 
mass. In other words, random velocities of molecular clouds of 
different masses are not expected to be in equipartition of en¬ 
ergy (Mdv^/2 =const.). Observations by Stark & Lee (2005) 
have shown that the random velocities of low-mass molecular 
clouds (< 2 X IO^MG) only vary by a few, with no dependence 
on cloud mass. These observations are therefore more consistent 
with our picture than a model in which molecular clouds acquire 
their relative velocities via mutual gravitational interaction. 

In limited circumstances, created molecular clouds collide 
with one another. This produces highly gravitationally unstable 
molecular gas an d may trigger very active star formation (e.g., 
iFukui et al.ll2014t) . Inoue & Fukui (2013) have done magnetohy- 
drodynamical simulations of a cloud-cloud collision and argue 
that it may lead to active formation of massive stars (see also 
IVaidva et al.ll2013h . This mode of massive star formation is not, 
however, a prerequisite of our model. 

2.1. Formation Timescale of Molecular Clouds 

Let’s first mo del _ the growth of molecular clouds. 

Ilnoue & Inutsuka! (l2012h have shown that we need multi¬ 
ple episodes of compression of HI clouds to create molecular 
clouds. According to the standard picture of supernova- 
regulated ISM dynamics (e.g., McKee & ©striker 1977), 
the typical timescale between consecutive compressions by 
supernova remnants is about IMyr. The total creation rate of ex¬ 
panding bubbles is larger than the occurrence rate of supernova 
explosions, since the former can also be created by radiation 
from massive stars less massive than supernova progenitors. 
Therefore, the actual timescale of compressions in ISM, Texp, 
should be somewhat smaller than IMyr if it is averaged over 
the Galactic thin disk. Obviously the compression timescale is 
smaller in the spiral arms and larger in inter-arm regions since 
star formation activity is concentrated in the spiral arms. Thus, 
we have to consider the time evolution of cloud mass for much 
longer than 1 Myr. 

Let us estimate the typical timescale of molecular cloud 
growth. Ilnoue & Inutsukal d2012h have shown that the angle be¬ 
tween the converging flow direction and the average direction of 
the magnetic held should be less than a certain angle for molec¬ 
ular cloud formation. Although Inoue & Inutsuka (2009) shows 
that this critical angle depends on the flow speed, we adopt a 
critical angle of 15 degrees (=0.26 radian) in the following dis¬ 
cussion for simplicity. This value is not so different from the an¬ 
gle (~ 20 degrees) for possible compression in the si mpler one¬ 
dimensional model by iHennebelle & Perauld (l2000li . For sim¬ 
plicity we assume that magnetic held is uniform in the region we 
consider and the direction of compression is isotropic. The solid 
angle spanned by the possible directions of compression result¬ 
ing in the formation of a molecular cloud is 0.26^71. The anti¬ 
parallel directions are also possible. Therefore, the probability, 
p, of successfully forming a molecular cloud in a single com¬ 
pression can be estimated by the ratio of solid angle over which 
compressions lead to molecular cloud formation to the solid an¬ 
gle of the whole sphere, i.e., p — 2 ■ 0.26^71/(4;7r) = 0.034. Note 
that Figure 1 is not the snapshot just after the birth of molec¬ 
ular clouds, but the result of one additional compression of the 
molecular clouds in which the direction of compression is per¬ 


pendicular to the mean direction of the magnetic held lines. We 
also emphasize that since the formation of a GMC requires many 
episodes of compression, our model does not predict a strong 
correlation between the present-day magnetic held direction and 
the orientation of the GMC. 

After each compression a cloud may slightly expand because 
of the reduced pressure of the ambient medium, which may re¬ 
sult in the loss of diffuse components of cloud mass. Observa- 
tionally the average column densities of molecular clouds do not 
seem to change very much and always appear to correspond to a 
visual extinction of several. This means that the mass of a cloud 
is proportional to its cross-section. Since the compressional for¬ 
mation of molecular material is expected to be proportional to 
the cross-section of the pre-existing cloud, we can model the rate 
of increase of molecular cloud mass as 

dM _ M 

~ ^ ’ 

where T{ denotes the formation timescale. This equation shows 
that resultant mass of each molecular cloud grows exponentially 
with a long timescale Tf if we average in time over a few Myr. 
If self-gravity increases the accumulation rate of mass into the 
molecular cloud, the right-hand side of Equation (1) may have a 
stronger dependence on mass. For example, the so-called “grav¬ 
itational focusing factor” increases the cross section of coales¬ 
cence by a factor proportional to the square of mass for the large 
mass limit. This will produce a significantly steeper slope of the 
cloud mass function (see Section 4). A linear dependence on 
mass in our formulation implicitly assumes that self-gravity of 
the whole molecular cloud does not significantly affect the cloud 
growth. 

Based on our investigation of molecular cloud formation de¬ 
scribed above, we estimate the formation timescale as follows: 

Ti - — ■ Texp. (2) 

P 

The average value in spiral arm regions would be Tf ~ 10 Myr, 
but can be factor of a few longer in the inter-arm regions. In re¬ 
ality, many repeated compressions with large angles between the 
flow direction and mean magnetic field lines gradually increase 
the dense parts of clouds, and hence contribute to the formation 
of molecular clouds over a long timescale. This may mean that 
the actual value of Tf is somewhat smaller than the estimate of 
Equation (2). 

Eukui et al. (2009) has shown that the clouds with masses 
(a few xIO^Mq) gain their mass at a rate 0.05 Mo/yr over a 
timescale lOMyr. This means that the mass of a cloud in their 
sample doubles in ~ lOMyr, which is consistent with our choice 
of Tf - lOMyr. Note, however, that Eukui et al. (2009) argued 
that the atomic gas accretion is driven by the self-gravity of a 
GMC, which is not included in the present modelling where we 
assume that gas accretion is essentially driven by the interac¬ 
tion with expanding bubbles. If the gravitational force is signif¬ 
icant for the HI accretion onto GMC, it possibly enhances the 
growth rate of molecular cloud (i.e., smaller Tf). Further quanti¬ 
tative studies of the effect of self-gravity on the accretion of gas 
onto a GMC remain to be done. In the present paper we neglect 
this effect and do not distinguish self-gravitationally bound and 
pressure-confined clouds, for simplicity. 

In regions where the number density of molecular clouds is 
very large, cloud-cloud collision may contribute to the increase 
of cloud mass, and hence, may also affect the mass function of 
molecular clouds. The detailed modelling of cloud-cloud colli¬ 
sion will be given in our forthcoming paper. Here we ignore the 
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contribution of cloud-cloud collision to the change of mass func¬ 
tion and simply use the constant value of Tf. 


3. Quenching of Star Formation in Moiecular 
Ciouds 

Next we consider the destruction of molecular clouds to de¬ 
termine how the star formation is quenched. Dale et al. (2012, 
2013) have done extensive three dimensional simulations of star 
cluster formation with ionization or stellar wind feedback and 
shown that the effects of photo-ionization and stellar winds are 
limited in quenc hing the star forma ti on in massive molecular 
clouds (see also IWalch et al. I l 2012 h . iDiaz-Miller et al.l (Il998h 
calculated the steady-state structures of HII regions and pointed 
out that the photodissociation of hydrogen molecules due to 
FUV photons is much more important than photoionization 
due to UV photons for the destruc ti on of molecular clouds. 
[Hosokawa & Inutsukal (l2005l l2006allbl l2007h actually included 
photodissociation in the detailed radiation hydrodynamical cal¬ 
culations of an expanding HII region in a non-magnetized ISM 
(by resolving photodissociative line radiation), and found the 
limited effect of ionizing radiation and essentially confirmed the 
critical importance of FUV radiation for the ambient molecular 
cloud. 



10 100 


stellar mass: M. [ Mg ] 


3.1. Expanding HII Regions in Magnetized ISM 


In the case of non-magnetized molecular gas of density lO^cm 


around a massive star larger than ~ 20Mq, a large amount of 
gas (~ 3 X 10 "^Mq) is photodissociated and re-processed into 
molecular material in the dense shell around the ex panding HII 


region within 5 Mvrs dHosokawa & Inutsukal2006b 

). According 

to the series of papers bv|lnoue & Inutsukal (l2008l 

20091). how- 


ever, the inclusion of magnetic field is expected to reduce the 
density of the swept-up shell substantially. Therefore the mag¬ 
netic field should affect significantly the actual structure of com¬ 
pressed shell and subsequent st ar formation process (c.f., 3D 
simulation bv I Arthur et ai.ll 20 l'lll . 


To quantitatively analyze the consequence, we have done 
numerical magnetohydrodynamics simulations of an expand¬ 
ing bubble due to UV and FUV photons from the massive 
star. The details of the meth od is the same as described in 
[Hosokawa & Inutsukal (l2006allbl) except for the inclusion of the 
magnetic field. Since the calculation assumes spherical symme¬ 
try, we include only the 13/rGauss magnetic field that is trans¬ 
verse to the radial direction as a simplification. The magnetic 
pressure due to transverse field is acc ounted for in the Riemann 
solver as in Sanoetal.1 ( 1999h (see ISuzuki & Inutsukal |2005 
llwasaki & Inutsuka 11201 1). 


The upper panel of Figure [3 shows the resultant masses of 
ionized gas in HII region and atomic gas in photo-dissociation 
region transformed from cold molecular gas around an expand¬ 
ing HII region at the termination time as a function of the mass 
of a central star (M*). Also plotted is the warm molecular gas 
in and outside the compressed shell around the HII region. The 
temperature of the warm molecular gas exceeds 50K. Its column 
density is smaller than 10 ^'cm^^, and hence, dust shielding for 
CO molecule is not effective and all the CO molecules are photo- 
dissociated. This warm molecular gas without CO (so-called, 
CO-dark H 2 ) is not expected to be gravitationally bound unless 
the mass of the parental molecular cloud is exceptionally large. 
Therefore the subsequent star formation in this warm molecular 


Fig. 3. Upper Panel: Masses in various phases transformed from cold 
molecular gas around an expanding HII region at the termination time 
as a function of the mass of a central star. The red dot-dashed line, the 
blue dotted line, and purple dashed line correspond to ionized hydro¬ 
gen in the HII region, neutral hydrogen in the photodissociation region, 
and warm molecular hydrogen gas without CO, respectively. The upper¬ 
most black solid line denotes the total mass of these non-star-forming 
gases. Lower Panel: The IMF-weighted mass of non-star-forming gas 
transformed from molecular gas by a massive star of mass M,. The area 
under the curve is proportional to the mass generated by massive stars 
whose mass function follows dnldflogM,) oc The peak of the 

curve determines the inverse of star formation efficiency esp (see expla¬ 
nation below Equation|7]l. 


gas is not expected. The uppermost black solid line denotes the 
total mass (Mg(M»)) of these non-star-forming gases. 

The lower panel of Figure 3 shows gas mass in the upper 
panel multiplied by ^ (blue dashed curve), (black 

solid curve), and ^ (red dotted curve). The areas under these 
curves are proportional to the mass affected by massive stars 
whose mass distribution follows dritldilogM,) oc 
and Ml’’, respectively. We can see that the shape of the curve 
does not vary much, and stars with mass 20 ~ 3QMq always 
dominate the disruption of the molecular cloud. 

Our calculations include ionization, photodissociation, and 
magnetohydrodynamical shock waves, but are restricted to a 
spherical geometry. Therefore we should investigate the dis¬ 
persal process in more realistic three dimensional simulations. 
However, the inclusion of photo-dissociation requires the nu¬ 
merical calculation of FUV line transfer and hence remains ex¬ 
tremely difficult in multi-dimensional simulations. 

3.2. Star Formation Efficiency 

Hereafter we assume the power law exponent of mass in the 
initial mass function of stars for large mass (M, > 8 M 0 ) is 
-yS and 2 < < 3 (c/wt/c/M* oc M,^). Now we calculate 
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the total mass of non-star forming gas disrupted by new born 
stars in a cloud. One might think that the total mass of non- 
star-forming mass in the stellar system can be calculated by 

nOO 

-^g,total = Jq M^(M,)(dn^ldM,)dMt. However, this estimation 
is meaningful only in the case the number of massive stars in a 
cloud are very large, f (dn^jdMt)dMt » 1. In reality, the 
number of massive stars in a molecular cloud of intermediate 
mass is quite small, and even a single massive star can destruct 
the whole parental molecular cloud. Thus, to analyze the quench¬ 
ing of star formation in molecular clouds, it is more appropriate 
to determine the most likely mass of the star that is responsible 
for the destruction of molecular clouds. 

Since we assume the large mass side of the stellar initial mass 
function can be approximated by the power law of the exponent 
-yS, we can express the mass function in logarithmic mass of 
massive stars created in a cloud as 


drit dn, / M, \ 

-= M,- = Nt - 

dlogM, dM, [MqI 


for M. > 8Mo 


( 3 ) 


Note that the pre-factor is defined for the mass distribution of 
stars in the individual cloud we are analyzing. For convenience, 
we define the effective minimum mass (M.m) of a star in the 
hypothetical power law mass function by the following formula 
for the total mass in the cloud. 



Suppose that a single massive star more massive than M*d is 
created in the molecular cloud. This condition can be expressed 
as 


/^OO 

•JM,i 


drif 

Im, 


■dM, = 


N, 

P-I 




B-I 


for M* > 8Mo. 


( 5 ) 


This equation relates M*d and N,. We can express the total mass 
of stars in the cloud as a function of M,d by eliminating Nt in 
equations 0 and 0, 


total 


r’ 

,/3-2j[M,^j [MoI 


for M. > 8Mo. 


(6) 


Thus, Mj totai ^ ^ 8Mq. Now we suppose that 

a molecular cloud of mass is eventually destroyed by UV 

and FUV photons from a star of mass M,d born in the cloud, and 
hence, star formation in the cloud is quenched. The condition for 
this to occur can be written as = Mg and esF^d = Af*,total, 
where esF is the star formation efficiency (the ratio of the total 
mass of stars to the mass of the parental cloud). 

If esF is smaller than the value that would satisfy the above 
condition, the cloud destruction is not sufficient and star forma¬ 
tion continues using the remaining cold molecular material in 
the cloud, which in turn increases esF- Thus, we expect that the 
actual evolution of a molecular cloud finally satisifies the above 
condition when the star formation is eventually quenched. This 
means that the star formation efficiency should be given by 

^ M,,total ^ 

Mg(M.d) \li-2)\M,J 

The preceding argument suggests that the value of star formation 
efficiency should take the minimum value of the right hand side 


of this equation, i.e., the maximum value of where Mg 

is a function of M,. Figure 3 shows that the maximum value of 
is attained at ~ 30Mo where Mg is about IO^Mq. 
Therefore we can conclude that once a massive star of M^d = 
30 + lOM© is created, the star formation is eventually quenched 
in a cloud of mass Mg(M,d) ~ 10^M©. This corresponds to esF ~ 
10“^, if we adopt fi - 2.5 and = O.IM©. The dependence 
of esF on Mtm is quite weak (/3 - 2 ~ 0.5) as shown in Equation 
(7). It is also not sensitive on JS in the limited range (2.3 < jd < 
2.7) as shown in Figure 3. Thus, the authors think that this value 
of esF ~ 10“^ is robust in typical star forming regions in our 
Galaxy. This argument may explain the reason for the low star 
formation efficiency in molecular clouds observati onally found 
many decades ago (e.g. JZuckerman & Evansiri974t) . 

Note that a sharp increase of is due to the sharp 

increase of UV/EUV luminosity at M, ~ 2QMq. Therefore a star 
much smaller than 20Mq is not expected to be the main disrupter 
of the molecular cloud. Eor example, the upper panel of Eigure 3 
shows that a lOM© star can quench ~ lO^M© of the surrounding 
molecular material. However a lO^M© molecular cloud is not 
likely to produce a lOM© star unless esp ~ 1 as can be seen in 
equation (O, and hence, the destruction of lO^M© cloud by a 
lOM© star is not expected, in general. 

If the initial mass function does not depend on the parent 
cloud mass as we assume here, the star formation efficiency is 
not expected to depend on mass for a cloud larger than ~ 10^M©. 
This can be understood as follows. The number of UV/EUV pho¬ 
tons is proportional to the number of massive stars, which in¬ 
creases with the mass of the cloud. However, the required num¬ 
ber of photons also increases with the mass of the cloud. Eor ex¬ 
ample, a lO^M© cloud will produce 10 stars with mass > 30M© 
if esF = 10^^, p — 2.5, and M.m = O.IM©. Then, these 10 
massive stars will destroy 10 x 10^M© molecular gas. Therefore 
star formation in the whole molecular cloud is quenched when 
esF = 10“^. Therefore we can conclude that the star formation 
efficiency does not depend on the mass of the cloud if the shape 
of the initial mass function does not depend on the mass of the 
cloud. 

Now we can estimate the timescale for the destruction of 
a molecular cloud. Our calculation of the expanding ioniza¬ 
tion/dissociation front in the magnetized molecular cloud shows 
that a ~ 10^M© molecular cloud can be destroyed within 4 Myrs. 
The actual destruction timescale, T^, should be the sum of the 
timescale of formation of a massive star and expansion timescale 
of the HII region, i.e., Td ~ r* -F 4Myr, where T, denotes the 
timescale for a massive star to form once the cloud is created. 
After one cycle of molecular cloud destruction over a timescale 
Td, only a fraction, esp, of the molecular gas is transformed 
into stars. Therefore, the timescale to completely transform a 
molecular cloud to stars is Td/esF ~ 1.4Gyr for T, ~ lOMyrs 
and Ti ~ 14Myrs. This may explain the so-called “depletion 
timescale” of molecular clouds th at manifests observationally in 
the Schmidt-Kennicutt Law (e. g., iBigiel et'^l201 ll : iLada et aTI 
l2012l : lKennicutt & Evansll2012h . 


4. Mass Function of Molecular Clouds 

In order to describe the time evolution of the mass function of 
molecular clouds, na{M) - dN^ldM, over a timescale much 
longer than 1 Myr, we adopt coarse-graining of short-timescale 
growth and destruction of clouds, and describe the continuity 
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equation of molecular clouds in mass space as 


dtici d i dM\ tici 


(8) 


where n^\{dMldt) denotes the flux of mass function in mass 
space, dMjdt describes the growth rate of the molecular cloud 
as given in Equation (1). The sink term on the right hand side 
of this equation corresponds to the destruction rate of molecular 
clouds in the sense of ensemble average. If the dynamical effects 
such as s hear and tidal stre s ses contribute to the clo ud destruc¬ 
tion te.g.. lKoda et al.ll2()()^ iDobbs & Pringlell2013h . we should 
modify Td in this equation. Since the left hand side of this equa¬ 
tion should be regarded as the ensemble average, the term l/Td 
represents the sum of the destruction rate of all the possible pro¬ 
cesses. Here we simply assume that the resultant Td is not very 
different from our estimate of destruction due to radiation feed¬ 
back from massive stars. 

According to the series of our work on the formation of 
molecular clouds, the molecular cloud as a whole is not neces¬ 
sarily created as a self-gravitationally bound object. Therefore, 
our modelling of the mass function of molecular clouds is not 
restricted to the self-gravitationally bound clouds. However, our 
modelling is not intended to describe the spatially extended dif¬ 
fuse molecular clouds much larger than the typical size of the 
bubbles (< lOOpc). 

A steady state solution of the above equation is 


n,,(M) = — 
Mo 



( 9 ) 


where Nq is a constant and 


1 

a = 1 H-. 

Ti 


( 10 ) 


For conditions typical of spiral arm regions in our Galaxy, we 
expect r* ~ Tf and thus Tf % Tj, which corresponds to 
1 < a < 2. For example, Tf - = lOMyrs corresponds to 

a » 1.7, which agrees well with observat i ons ([Solomon et al 


j_987: [&amer et akll998t[Hever et afl200 it [Roman Duval et al 


20101 ) . 

However, in a quiescent region away from spiral arms or in 
the outer disk, in which there is a very limited amount of dense 
material, Tf is expected to be larger at least by a factor of a few 
than in spiral arms. In contrast, Td is not necessarily expected to 
be large even in such an environment, since the meaning of Td 
is the average timescale of cloud destruction that occurs after 
the cloud is created, and thus, it does not necessarily depend 
on the growth timescale of the cloud. Therefore, we expect that 
Td can be smaller than Tf in such an environment, which may 
produce a - 1 + Tf/T^ > 2. This tendency is actu ally observed 
in Milky Way outer disk, FMC, M33, and M51 dRo solowskvl 


2005tlWong et al.ll201lUGratier et al.ll2^lHughes et al.ll2010 : 

Colombo et al.ll2014l) . 

The total number of molecular clouds is calculated as 


^total 


nM: 

JMi 


n{M)dM = 


No 


/Mo 

a - 1 

{MfJ 

[m2 


No /Mo 


a - 1 \ Ml 


( 11 ) 


where we used M 2 » Mi in the final estimate. The total number 
of clouds is essentially determined by the lower limit of the mass 


of the cloud. Fikewise, the total mass of the molecular clouds is 

\ 2~a' 

I l\f\IVI^ I IVI'} \ I IV /1 


/ 


Ml 


a 


»-c \ 2—0' 

M2\ 

/Ml 

Mo/ 

Uo, 


AoMp /MzY 

2-a \Mo/ 


( 12 ) 


where we used M 2 » Mi and a < 2 in the final estimate. Thus, 
the total mass of molecular clouds is essentially determined by 
the upper limit of the mass of the cloud. Fet us assume Mtotai ~ 
10 ^Mo in the Galaxy, then our simple choice of Mi = lO^M©, 
M 2 = 10®Mq, and a = 1.5 corresponds to Atotai ~ 10® and 
the ayerage mass of molecular clouds is Mave = IHtotai/AYotai ~ 
IO^^Mq. Note that these numbers depend on our choice of Mi. 


5. Summary 

In general, dense molecular clouds cannot be created in shock 
wayes propagating in magnetized WNM without cold HI clouds. 
In this paper we identify repeated interactions of shock wayes 
in dense ISM as a basic mechanism for creating fllametary 
mole cular clouds, which are ubiquitously obseryed in the nearby 
ISM (I Andre et al.l 12014) . This suggests an expanding-bubble- 
dominated picture of the formation of molecular clouds in our 
Galaxy, which enables us to enyision an oyerall picture of the 
statistics of molecular clouds and resultant star formation. To¬ 
gether with the findings of our preyious work, our conclusions 
are summarized as follows; 

1. Turbulent cold HI clouds embedded in WNM can be read¬ 
ily created in the expanding shells of HII regions or in the 
yery late phase of supernoya remnants. In contrast, the for¬ 
mation of molecular clouds in a magnetized ISM needs many 
compression eyents. Once low-mass molecular clouds are 
formed, an additional compression creates many filamentary 
molecular clouds. One compression corresponds to of order 
IMyr on ayerage in our Galaxy. The timescale of cloud for¬ 
mation is a few times lOMyrs. 

2. Since the galactic thin disk is occupied by many bubbles, 
molecular clouds are formed in the oyerlapping regions of 
(old and new) bubbles. Howeyer, since the ayerage lifetime 
of each bubble is shorter than the timescale of cloud for¬ 
mation, it is difficult to obseryationally identify the multiple 
bubbles that created the molecular clouds. 

3. The yelocity dispersion of molecular clouds should originate 
in the expansion velocities of bubbles. This is estimated to 
be <10km/s and should not strongly depend on the mass of 
the molecular cloud. 

4. To describe the growth of molecular cloud mass we can tem¬ 
porally smooth out the evolution over timescales larger than 
~ IMyr. The resultant mass (smoothed over time) of each 
molecular cloud is an almost exponentially increasing func¬ 
tion of time. 

5. The destruction of a molecular cloud is mainly due to 
UV/FUV radiation from massive stars more massive than 
2QMq. The probability of cloud destruction is not a sensi¬ 
tive function of the mass of molecular clouds. If the shape of 
the initial mass function does not vary much with the mass 
of parent molecular clouds, cloud destruction by 30 + IOM 0 
stars results in a star formation efficiency of order 1%. This 
property explains the observed constancy of the gas deple¬ 
tion timescale (1 ~ 2 Gyr) of giant molecular clouds in the 
solar neighborhood and possibly in some external galaxies 
where the normalizations for the Schmidt-Kennicutt Faw ob¬ 
tained by high-density tracers are shown to be similar. 
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6 . The steady state of the evolution of the cloud mass function 
corresponds to a power law with exponent -n in the range 
1 < n < 2 in the spiral arm regions of our Galaxy. However, 
a larger value of the exponent, such as n > 2, is possible in 
the inter-arm regions. 

Note that the first and th ird conclusions have par tly shown in 
our previous investigations (llnoue & lnutsukal2009h . In addition 
we can suggest the following implications from these conclu¬ 
sions: 

7. Star formation starts, even in small molecular clouds, once 
the line-mass of a n internal self-grav itating filament exceeds 
the critical value (lAndre et alj|2010l) . Our analysis suggests 
that the mass of an individual molecular cloud increases 
roughly exponentially over ~ 10 Myrs. According to the for¬ 
mation mechanism driven by repeated compressions, we ex¬ 
pect that the total mass in filaments of sufficiently high line- 
mass increases with the number of compressional events. 
This means that the mass of star-forming dense gas increases 
with the mass of the molecular cloud and the star formation 
should accelerate over many million years. This conjecture 
may provide a c lue in understanding th e star formation his¬ 
tories found by IPalla & Stahled (l2000l) in seven individual 
molecular clouds such as Taurus-Auriga and p Ophiuchi. 

8 . Molecular clouds may collide over a timescale of a few times 
10 Myrs, depending on the relative locations in adjacent (al¬ 
most invisible) bubbles. Such a molecular cloud collision 
may result in active star formation in a giant molecular cloud 
Ic.g.- lpukui et ^l2014l) . 

These implications should be investigated in more detail 
by numerical simulations. Our radiation magnetohydrodynamics 
simulations of an expanding bubble due to UV and FUV photons 
from the massive star show that most of the material in molecular 
clouds become warm molecular clouds without CO molecules. 
Although we have to investigate the fate of the CO-dark gas in 
more detail, we expect that the total mass can be very large and 
may account for the dark gas indicated by various observations 
(e.g.. lGrenier et al.ll2005h . 

There are many report that the Kennicutt-Schmidt cor- 
relation varies with some properties of ga l axies (e.g., 
Saintonge et al.l 1201 iL 120121: iMeidt et ^ 1201 3t iDavis et ^ 
2014h . In addition , a simple relation does not fit to th e center of 
our Galaxy (e.g.. iLongmore. Ballv. Testj et al.|l2013h . The rea¬ 
sons for these deviations remain to be studied. 

Acknowledgements. SI thanks Hiroshi Kobayashi and Jennifer M. Stone for use¬ 
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